Hydrogen-induced ferromagnetism in ZnO single crystals investigated by 

Magnetotransport 
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We investigated the electrical and magnetic properties ol low-energy hydrogen-implanted ZnO 
single crystals with hydrogen concentrations up to ~ 3 at. % in the first 20 nm surface layer between 
10 K and 300 K. All samples showed clear ferromagnetic hysteresis loops at 300 K with a saturation 
magnetization up to ~ 4 emu/g. The measured anomalous Hall effect agrees with the hysteresis 
loops measured by superconducting quantum interferometer device magnetometry. All the H-treated 
ZnO crystals exhibited a negative magnetoresistance up to the room temperature. The relative 
magnitude of the anisotropic magnetoresistance reaches 0.4 % at 250 K and 2 % at 10 K, exhibiting 
an anomalous, non-monotonous behavior and a change of sign below 100 K. All the experimental 
data indicate that hydrogen atoms alone in a few percent range trigger a magnetic order in a ZnO 
crystalline state. Hydrogen implantation turns out to be a simpler and effective method to generate 
a magnetic order in ZnO, which provides interesting possibilities for future applications due to the 
strong reduction of the electrical resistance. 

PACS numbers: 75.30.Cr, 75.50.Pp, 75.60.Ej 



I. INTRODUCTION 

Defect-induced magnetism (DIM) appears to be a gen- 
eral phenomenon observed in nominally non-magnetic 
solids starting from the archetype graphite [11, to sev- 
eral oxides like ZnO, pure or doped with non-magnetic 
elementsM, HfOz I, TiOz E El, SrTiOg Bill, 
SrO:N [iSj as well as Si-based samples [T^. to mention 
only a few examples from a large number being reported 
nowadays (for recent reviews on this subject see Refs.lTsl- 
ITsI ). Experimental facts demonstrate that defects, like 
vacancies, without or with the presence of non-magnetic 
ad-atoms, play a main role in triggering magnetic order 
in these systems. Recently, room temperature ferromag- 
netism was reported in Cu-doped ZnO films, investigated 
by soft x-ray magnetic circular dichroism [l9[. The re- 
sults of this study strengthens the existence of the DIM 
phenomenon in general and in ZnO in particular. 

ZnO is a wide band gap semiconductor, which crystal- 
lizes in the hexagonal wurtzite, zincblcnde and rocksalt 
structures. However, hexagonal wurtzite is the most in- 
tensively studied crystal structure of ZnO because of its 
potential applications in the field of spintronics, trans- 
parent electronics, piezoelectricity, optoelectronics, etc. 
Hydrogen is one of the most abundant and unavoidable 
impurities in ZnO. The presence of hydrogen can infiu- 
ence the electrical and the magnetic properties of ZnO. 
The role of hydrogen in enhancing the ferromagnetism in 
3d-transition metal-doped ZnO was recently studied ex- 
perimentally [l^l and theoretically dlj. On the other side 
there are theoretical studies reporting on the possibility 
of room temperature ferromagnetism due to hydrogen 
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adsorption at the surface of ZnO [2^ |23| . We note how- 
ever that there is a lack of experimental and theoretical 
studies on the possibility of hydrogen-induced magnetic 
order inside the crystalline structure of pure ZnO. 

In this work we are interested on DIM in ZnO crys- 
tals triggered through the implantation of hydrogen at 
low energies and its influence on the magnetotransport. 
We report on a detailed experimental study that demon- 
strates how the intentional doping of hydrogen in the 
percent range and at low-enough energies influences sub- 
stantially the electrical, magnetic and magnetotransport 
properties of ZnO single crystals. In particular, we show 
in this report a hydrogen-induced anisotropic magne- 
toresistance (AMR) as well as the anomalous Hall effect 
(AHE) in H-ZnO single crystals. Therefore, the existence 
of hydrogen-induced ferromagnetism in H-ZnO samples 
is supported not only by the usual magnetization data 
taken with a superconducting quantum interference de- 
vice (SOUID)[24| but also from magnetotransport mea- 
surements. 



II. EXPERIMENTAL DETAILS 

Hydrothermally grown ZnO (0001) single crystals of di- 
mensions (6x6x0.5) mm'^ with both sides polished were 
supplied by CrysTec GmbH. The ZnO samples were ex- 
posed to remote hydrogen dc plasma for different time in- 
tervals in a parallel-plate system. The voltage difference 
between the two plates was kept at 1 kV. The samples 
were mounted on a heater block held at a fixed temper- 
ature of 400 °C and they were placed ~ 100 mm down 
stream from the plasma with a bias voltage of ~ -330 V. A 
bias current of ~ 50 fiA was measured during the plasma 
treatment. The pressure in the chamber during the pro- 
cess was maintained around 1 mbar. Three samples H-1, 
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H-2 and H-3 for time intervals 30, 60 and 90 min, re- 
spectively, were treated in the H-plasma chamber. The 
implantation depth (for the chosen energy) as well as for 
the concentration characterization analysis (see below) 
were estimated using SRIM[25|. From this Monte Carlo 
simulation program we estimate a penetration depth of 
20 nm for the implanted hydrogen atoms. As experi- 
mentally shown in Rcf . [2J] the main ferromagnetic signal 
comes from this near surface region, in agreement with 
the estimates. 

Nuclear reaction analysis (NRA) was used to deter- 
mine the hydrogen concentration in ZnO crystals before 
and after hydrogen plasma treatment [2^. The NRA 
has a depth resolution of 5 nm with an average error 
in the concentration of 0.02%. The hydrogen concentra- 
tion in ZnO crystals measured by NRA before and after 
remote hydrogen treatment was found to be 0.14±0.03 
and 0.64±0.07 at.% in the first 200 nm from the sur- 
face, respectively [l^. From this concentration analysis 
we conclude that the first 20 nm near surface region 
should have a hydrogen concentration of the order of 
~ 3 at.% for sample H-3. The hydrogen concentrations 
for samples H-1 and H-2 are ^ 1 at.% and 1.8 at.%, 
respectively. 

We performed Particle Induced X-ray Emission 
(FIXE) measurements to analyze the concentration of 
magnetic impurities in the H-ZnO samples. There was no 
significant difference in the Fe concentration (< 60 ppm) 
before and after H-plasma treatment. 

The magnetization of the ZnO single crystals before 
and after the plasma treatment was determined with a 
SQUID. The magnetotransport measurements were per- 
formed in an Oxford cryostat with a magnetic field up to 
8 T and a rotating sample holder allowing us to measure 
the resistance at different angles between magnetic field 
and the input current. The electrical contacts on the 
samples were prepared using silver paste in a Van der 
Pauw configuration. The I /V characteristics were mea- 
sured to check for deviations from the ohmic behavior. 
All the transport data presented in this paper were taken 
in a linear, ohmic regime. The resistance was measured 
with an ac resistance bridge with a relative resolution of 
0.01%. 



III. RESULTS AND DISCUSSION 



and the measured signal is given by the addition of the 
diamagnetic plus the ferromagnetic one of the near sur- 
face region. After subtraction of the diamagnetic con- 
tribution the magnetization coming from the ferromag- 
netic part was calculated assuming a ferromagnetic mass 
homo gene ously located at the first 20 nm near surface 
rcgion[2J|, see Fig. [Ijb). We see that the magnetization 
increases with H-implantation as reported in Ref. [l^ and 
also decreasing temperature, compare the results of sam- 
ple H-3 in Fig. [TJb) with those in the inset of Fig. [Ija). 
The ferromagnetic magnetization per total volume of a 
virgin, untreated ZnO crystal of the same type as used 
here, is of the order of 10~^ emu/gQ. For comparison 
and taking into account that a similar near surface re- 
gion in nonmagnetic oxide crystals could be the source 
for the ferromagnetic signals Q, the untreated ZnO crys- 
tal would have a saturation ferromagnetic magnetization 
half of that of sample H-1. 

All samples exhibit a coercivity of 18-20 mT at room 
temperature. The saturation magnetization increases by 
increasing the H-concentration and reaches ~ 4 emu/g 
for sample H-3. The ferromagnetic magnetization at 5 K 
increases by ~ 50% of its value at room temperature 
(see inset in Fig. 1(a)). From the measurement of the 
ferromagnetic remanent magnetic moment we estimate a 
Curie temperature of 450 ± 25 K [24 1. 

Recently, N. Sanchez et al.[2^ have theoretically 
shown that atomic hydrogen adsorbed on the Zn-ZnO 
(0001) surface can form strong H-Zn bonds and lead to 
a metallic surface with a net magnetic moment of 0.5 fiB 
per H-atom. The obtained magnetization values indicate 
that our H-3 sample, for example, would have a net mag- 
netic moment of 0.2 per hydrogen atom. This esti- 
mate is obtained taking into account the amount of H im- 
planted. The difference between the two estimates may 
indicate that the assumed ferromagnetic mass is larger 
than the true one. This appears plausible because the 
hydrogen atoms are not necessarily homogeneously dis- 
tributed inside its penetration depth. We believe that in 
the first 20 nm depth we have regions where the magnetic 
order is less developed and therefore one tends to overes- 
timate the ferromagnetic mass. This picture of a rather 
inhomogeneous mixture of magnetic and non-magnetic 
regions is of importance to interpret the transport data, 
as we shall discuss below. 



A. Magnetization Measurements 

Figure [ija) shows the magnetic moment after subtrac- 
tion of the linear diamagnetic contribution vs. applied 
field measured at 5 K for all three samples with the 
SQUID. This signal is composed by two contributions. 
The main one is paramagnetic and follows the usual Bril- 
louin or Langevin function. If one subtracts it from the 
data a small ferromagnetic contribution still remains, as 
can be seen in the inset of Figure [TJa) for sample H-3. 

At 300 K the paramagnetic contribution is negligible 



B. Resistivity Measurements 

We measured the temperature dependence of the re- 
sistance of the three H-ZnO samples from 10 K to 
270 K, see Fig. [5] The resistance decreases increasing H- 
concentration, therefore it is reasonable to assume that 
the resistivity of the implanted part is much smaller than 
the resistivity from the rest of the single crystal, as has 
been already reported [28| . This assumption is supported 
by the direct comparison of the estimated values of the 
resistivity - assuming conduction within the 20 nm im- 
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FIG. 1: (a) Magnetic moment of the three H-ZnO crystals as 
a function of magnetic field at 5 K. The inset shows a ferro- 
magnetic contribution which is left after subtracting a para- 
magnetic one obtained by using the Langevin function, (b) 
Ferromagnetic magnetization of the same crystals at 300 K 
assuming a homogeneously distributed ferromagnetic mass at 
the first 20 nm near surface region [2^. A diamagnetic slope 
was subtracted from the data. Note that the magnetization 
of H-ZnO crystals increases with H-concentration. 



planted thickness, right y— axis in Fig. [5]- with the resis- 
tivity of the virgin ZnO single crystal, which is already at 
room temperature several orders of magnitude larger [29|. 

Before discussing in detail the observed behavior of 
the resistivity, we would like to describe shortly what is 
known about the hydrogen contribution to the formation 
and/or modification of the electronic band structure 
in ZnO. Hydrogen forms shallow donor states in bulk 
ZnO and is regarded as a source for n-typc conductivity. 
These shallow donor states are formed approximately 
30-60 meV below the conduction band. Upon doping 
level the donor energy states can be dispersed into an 
impurity band because of the Coulomb fields arising 
from the compensating acceptors and ionized donors [30j. 
This impurity band could further split into two bands, 
i.e. a lower band I and an upper band I~, which are 
formed with sirigle charged donors and neutral donors, 
respectively [Sll, 111] . 

As we have already noted in section 3A, we as- 
sume that the implanted H"*'-ions arc not homogeneously 
distributed in H-ZnO samples. Then, we may have 
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FIG. 2: Resistance of the three H-ZnO samples as a function 
of temperature at a zero applied field. At the right y— axis 
we show the resistivity estimated taken into account 20 nm 
implanted thickness of the single crystal. The inset blows out 
the high temperature part where the resistance shows a qual- 
itative change in its temperature dependence. The observed 
minimum shifts to the low temperatures by increasing the hy- 
drogen concentration. The red solid lines show the fits of the 
data to Eq. (2). The obtained values for Ri, R2 and R3 are 
(0.022, 0.014, 0.015) fl, (0.587, 0.612, 0.365) Q, and (197, 143, 
169) SI for samples (If-1, H-2, II-3), respectively. The rest 
resistance for all samples is Ro < lO^'' fl. 



hydrogen-rich metallic regions with resistance Rm{T) — 
Ro + RiT, embedded in a doped semiconducting matrix 
with resistance Rs{T) = R2Cxp{AE/2kBT) {ks is the 
Boltzmann constant). Rq, Ri,R2 arc free parameters as 
well as the activation energy AE which will be obtained 
by fitting the experimental data. Note that we consider a 
linear T dependence for the metallic part in all the tem- 
perature range. The metallic region contributes mainly 
at high enough temperatures, see inset in Fig. ^ and 
therefore it is unnecessary to assume a more complicated 
T— dependence that may be applicable at temperatures 
T < 100 K. We consider that these two contributions, 
the metallic- and semiconducting-like are in series. Due 
to the implantation distribution curve (2^ , it is clear that 
below ^ 20 nm a third intermediate region should exist 
that contributes with a resistance Rh (T) in parallel to the 
other two. The best fits have been achieved by assuming 
a variable range hopping (VRH)-likc mechanism as: 



Rh = -R3 cxp 



T 



1/5 



(1) 



Where R^ is a free parameter and i?„„ is a hopping en- 
ergy. The total resistance will be given then by 



i?(T)= {Rh{T)) ^ + (RrniT) + Rs{T)) ' 



(2) 



One can also take into account a fourth parallel re- 
sistance contribution arising from the pure ZnO single 



4 



crystal below a ^ 100 nm thick layer. However, the re- 
sistance of such pure ZnO crystal at all temperatures is 
in the range of Mfi and therefore its contribution to the 
total resistance is negligible. The fittings to the data of 
the three crystals are reasonably good as we can see in 
Fig. [2] The activation energy obtained from the fittings 
for the three H-ZnO single crystals is 60±2 meV. 

Qualitatively the observed temperature dependence of 
the resistivity is rather simple to understand. At temper- 
atures below 50 K, the resistance of the semiconducting 
contribution Rs is larger than the resistance Rh{T), be- 
coming this last the dominant transport mechanism [ssj 
because the thermal energy is not enough to excite the 
electrons from the upper impurity band I~ to the conduc- 
tion band. The hopping energy obtained from the fitting 
of the experimental data is Enn — 3 ± 0.5 meV. As the 
temperature increases the resistance of the H-ZnO sam- 
ples decreases following a semiconducting behavior with 
en effective activation energy AE ~ 60 mcV. The larger 
the hydrogen doping the lower is RmiT) and therefore 
the lower is the temperature of the minimum, see inset 
in Fig. [21 Although with this simple model we can un- 
derstand qualitatively the behavior measured in the re- 
sistance as a function of temperature and hydrogen con- 
centration, it does not provide us with a clear hint about 
the regions that contribute to the magnetic signal. If the 
magnetic order is confined mostly within the first 20 nm 
surface region[23| we expect then that cither the metallic 
or the semiconducting regions or even an intermediate 
region between these two and the VRH part, contributes 
to the ferromagnetic signal. As the semiconducting con- 
tribution overwhelms in a wide temperature range the 
metallic one, it appears plausible that this one may be 
responsible for the magnetoresistance behavior we have 
observed in H-ZnO single crystals. 



C. Charge Carriers 

We performed Hall measurements in a Van der Pauw 
configuration in order to obtain the charge carriers den- 
sity of our H-ZnO samples. Our Hall measurements con- 
firm that the conductivity in H-ZnO single crystals is n- 
type. We note that the estimated carrier concentration 
is an effective one obtained using the simplest expression 
n = I/Rhg, where Rh is the Hall resistance, assuming 
that only the sample volume of the first 20 nm H-rich 
surface layer of the ZnO single crystals contributes. In 
case electrons and holes would contribute with different 
densities and scattering rates then the two-band model 
is necessary to obtain the carrier densities. The use of 
its equations, however, implies the introduction of free, 
unknown parameters, like the scattering rates. In order 
to facilitate the comparison of our data with literature 
values we prefer to discuss the carrier density from the 
Hall data stressing that the measurable quantity is Rh- 
The temperature dependent of the carrier density tih is 
shown in Fig. |31 
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FIG. 3: (a) Carrier density of the three H-ZnO single crystals 
as a function of temperature. The carrier density shows an 
anomalous behavior around ~ 50 K. The solid lines are fit to 
the expression given in Eq. ( [S]) . 



The carrier concentration of samples H-1, H-2 and 



H-3 at room temperature are 6.72 xlO cm 



.56 



xlO^°cm ^ and 1.27 xlO^^cm ^. These values are com- 



parable to those found in, e.g., Ga-doped ZnO systeml34 1 
Above r ^ 50 K the carrier density increases with tem- 
perature following the equation: 



5cxp 



f AE \ 



\2kBT J ' 



(3) 



where a, b are free parameters. The activation energy ob- 
tained from the fits is AE = 60±2 meV, in agreement to 
the activation energy values obtained from the resistivity 
measurements, see Fig. O As expected, the carrier con- 
centration increases with hydrogen concentration. The 
increase in uh between the samples agrees roughly with 
the estimated increase in the hydrogen concentration. 

At temperatures T < 50 K, uh increases with a de- 
crease in temperature, see Fig. |31 This is an anomalous 
behavior that appears to be related to the change of the 
main contribution to the measured resistance, i.e. from 
the semiconducting region above 50 K to the VRH one 
below it, see Fig. [5J In this case it might be that the 
simple relation to estimate nniT < 50K) from the Hall 
resistance is not adequate and a more complicated equa- 
tion for the Hall signal of a material with two contribu- 
tions in parallel should be usedfssj. We note that the 
anomaly at T ^ 50 K is observed in all the magneto- 
transport properties we have measured, as we will show 
in the following sections. 



D. Magnetoresistance Measurements 

Apart from the magnetization measurements a further 
and important way to check for the existence of mag- 
netic order is through the measurement of the magne- 
totransport properties. Unlike magnetization measure- 
ments, magnetotransport properties are much less sensi- 
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FIG. 4: Magnetoresistance (in percentage) as a function of 
applied magnetic field at 10 K and 250 K for the three H-ZnO 
samples. A clear correlation between magnetoresistance and 
H-concentration is observed. The magnetoresistance increases 
with H-concentration. 
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FIG. 5: Magnetoresistance of sample H-1 (the other two sam- 
ples show a similar dependence) as a function of applied mag- 
netic field measured at several temperatures. The magnetore- 
sistance shows a negative temperature dependence but with 
an anomaly around 50 K. The solid (red) lines through the 
data points are fits to Eq. (Q. 



tivc to magnetic impurities, in case they remain below 
0.1%, and in general they reflect intrinsic characteris- 
tics of the sample. In this section we discuss the longi- 
tudinal magnetoresistance of H-ZnO where the magnetic 
field is applied parallel to the input current as well as to 
the sample main plane. The longitudinal magnetoresis- 
tance for the three H-ZnO samples measured up to 8 T 
at 10 K and 250 K is shown in Fig. |4l The magnetore- 
sistance is defined as [[R{H) - i?(0)]/i?(0)] where R{H) 
and i?(0) are the resistances with and without an applied 
magnetic field, respectively. 

All samples show a negative magnetoresistance at all 
temperatures and magnetic fields applied parallel to the 
main plane of the samples. The negative magnetoresis- 
tance has been observed in several other ZnO systems 
that show some kind of magnetic order (36ll37|. 

Figure [S] shows the temperature dependence of the 
magnetoresistance for the H-1 sample. Similar temper- 
ature dependent magnetoresistance is observed for the 
other two samples. It is clear from Fig. [5] that the mag- 
netoresistace of H-ZnO decreases in general with temper- 
ature. However, its field curvature at T ~ 50 K changes 
and from 50 K to 75 K the magnetoresistance increases 
with temperature. The decrease in the magnetoresis- 
tance with temperature is expected because the mag- 
netization at saturation of these samples also decreases 



with temperature, see Fig. [TJ As we have observed in 
the carrier concentration, the anomalous behavior in the 
magnetoresistance around 50 K might be related to the 
change of the main contribution to the resistance. Tak- 
ing into account that the main contribution to the mag- 
netic signal comes from the 20 nm surface contribution, 
the variable range hopping part and the magnetic semi- 
conducting part might have a common interface, which 
shows magnetic order and contributes to the magnetore- 
sistance at low temperatures. 

To elucidate our experimental results for the magne- 
toresistance we use a model proposed by Khosla and 
Fischer [s^ that combines negative and positive mag- 
netoresistances in semiconductors taking into account a 
third-order expansion of the s — d exchange Hamiltonian. 
The semiempirical formula is: 



= -a^ In 1 + IP-B^) + — - 

Po 1 + ^2^2 



(4) 



where c and d depend on the conductivity and the car- 
rier mobility, respectively. We will consider these as free 
parameters, and 



a' ^ AiJpF[S{S+l) + {M'^)l 



(5) 
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FIG. 6: Temperature dependence of the fitting parameters a, 
b, c and d of three H-ZnO samples. 
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(6) 



where /i is the mobility, g is the Lande f/-factor, a is 
a numerical constant. The fitting parameters a and b in 
Eq. (U) depend on several factors such as a spin scatter- 
ing amplitude Ai, the exchange integral J, the density 
of states at the Fermi energy pp, the spin of the local- 
ized magnetic moments S and the average magnetiza- 
tion square (M^). The negative first term in Eq. is 
attributed to a spin dependent scattering in third order 
s — d exchange Hamiltonian while the positive part (sec- 
ond term in the expression of Eq. (U)) takes into account 
field induced changes due to the two, s and d, conduction 
bands with different conductivities. 

The fits of the experimental data to Eq. (|4]) are shown 
in Fig. [5l The data can be well fitted with this model 
at all measured temperatures. All four fitting parame- 
ters show different temperature dependence (see Fig. 6). 
Note that the positive magnetoresistance in our samples 
is compensated by the large negative magnetoresistance 
contribution. Therefore, the uncertainty of the param- 
eters c and d is rather large because they are not in- 
dependent of the fitting procedure. Therefore, we con- 
centrate on the negative scattering contributions a and 
b. The fitting parameter b for the three H-ZnO sam- 
ples are plotted as a function of inverse temperature in 
Fig. [HI The parameter b shows a linear dependence in 
good agreement with theory. We found that the param- 
eter a is almost temperature independent and increases 
with H-concentration in the temperature range 50 K > 
T > 100 K indicating that the sample with higher H- 
concentration (H-3) is more magnetic, in agreement with 
the SQUID data shown in Fig. [TJ 

Parameters a and b defined in Eq. ([5]) and Eq. ^ re- 
spectively, are used in order to obtain the values of Jpp 
and Ai. The values of Jpp and Ai obtained from the 
experimental data are 0.56 and 0.14 for S = 1/2 and p 



= 36 cm^/V.s at 10 K, respectively. The values of Jpp 
and Ai for S = 3/2 are 0.33 and 0.12, respectively. The 
value of Jpp = 0.33 for S = 3/2 in H-ZnO is similar to 
the one obtained in CdS system Jpp = 0.4 [l^. These re- 
sults strongly suggest the contribution of s-d interaction 
in the H-ZnO system. 



E. Anisotropic Magnetoresistance 

There are two other magnetotransport effects, which 
are observed in our H-ZnO crystals and are worth men- 
tioning. One of them is the anisotropic magnetoresis- 
tance (AMR) effect. This effect represents the change in 
the resistance of a ferromagnetic material with the angle 
between the input current and applied field in plane. It 
is commonly associated with the presence of a spin split- 
ting of the electronic band at the Fermi level and a finite 
spin-orbit (L — S) coupling. The AMR arises in second 
order in the L — S coupling, in contrast to the magne- 
tocrystallinc anisotropy. In ferromagnetic materials with 
s— and rf— bands the AMR is understood arguing that 
the spin-orbit scattering increases the resistance by al- 
lowing a spin-flip and through this the occupation of free 
d— states in the corresponding spin dependent band. In 
general it is expected that the resistance is larger when 
the applied field is parallel to the current. We define the 
AMR amplitude as 



Ai?/i?avg — 



|fl(g||)-i?(g^)| 
iR{H^\) + R{H^))/2- 



(7) 



For polycrystalline ferromagnetic samples the change of 
the resistance due to the AMR has the following angle 
dependence: 



AR 



R{H = 0) 



= A cos 



(8) 



where 6 is the angle between the current / and the applied 
field or magnetization direction (in saturation) and A is 
a constant that depends on the density of d-states at the 
Fermi level, on the magnetization as well as on the sample 
quality. 

Figure [7] shows the field dependence of the magnetore- 
sistance for 0^0° and 90° at 10 K for the three H- 
ZnO samples. The angle 9 is the angle between the ap- 
plied current and the magnetic field. We observe a clear 
AMR effect in the three H-ZnO samples. We note that 
the AMR effect was already measured in ZnO but doped 
with Cojl^. Our results show clearly that a fundamen- 
tal property of ferromagnetic materials as the AMR can 
also be obtained by DIM in an oxide. The AMR effect 
in H-ZnO single crystals increases with H-concentration, 
see Fig. H 

However, for ferromagnetic single crystalline systems 
some influence from the lattice anisotropy can be ex- 
pected in the AMR and in this case the angular- 
dependent magnetoresistance is described by a Fourier 
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FIG. 7: Magnetoresistance of the H-ZnO samples at 10 K for 
two configurations of angles between current and the applied 
magnetic field. For clarity the data curves of the samples H-2 
and H[-3 are shifted by a constant value. 
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FIG. 8: Angular magnetoresistance of sample H-1 at a con- 
stant field of 5 T at several temperatures. The red solid lines 
are the fits to experimental data following Eq. 



series of cos{n6) and sm{n9) as 
AR 



R(H = 0) 



s„ sin(n6') + ^ c„ cos{n9) , (9) 



n=l 



where the coefficients s„ and c„ are related to the Hall 
and magnetoresistance contributions in the system. 



In order to study the AMR in more detail, we mea- 
sured the angular magnetoresistance at a constant field 
of 5 T at several temperatures, see Fig. [H] The variation 
of the angle between the magnetic field and the current 
ranges from —50° to 140°. This range of angle is nec- 
essary because we expect a 180° periodicity. There are 
several distinct features in the experimental data shown 
in Fig. [i (a) At T < 50 K the AMR is higher at 61 = 
° (field and current parallel to each other) than at 6 
= 90°, in agreement with the usual behavior, (b) At 
temperatures T > 50 K the AMR changes sign (or it 
shows a 90° shift in the angle dependence), see Fig. |S1 
(c) Figure [9] shows the AMR amplitude, see Eq.©, as 
a function of temperature for the three samples. There 
is a clear anomalous increase with temperature between 
50 K and 100 K, in spite of the fact the the magnetization 
at saturation decreases monotonously in all the temper- 
ature range according to the SQUID measurements (not 
shown). Note that the magnitude of the AMR in our 
H-ZnO samples at 250 K is ^ 0.4 %, a value compara- 
ble to the AMR observed in, e.g., Co films [ioj or Co:Cu 
multilayered nanowires[4l]|. (d) The measured angle de- 
pendence does not follow Eq. ^ applicable for polycrys- 
talline materials and indicates that the magnetic contri- 
bution in our H-ZnO single crystals comes from a single 
crystalline phase after H-plasma treatment. The AMR 
curves obtained for all three H-ZnO samples are fitted 
by Eq. ([9|) and the results of these fits are shown in 
Fig. |5] as (red) solid lines. The experimental data can 
be fitted quite well at all temperatures after expanding 
Eq. (jni) up to eighth order. The coefficients s„ and c„ are 
related to the antisymmetric (Hall) and symmetric (mag- 
netoresistance) contributions of the sample. The values 
of the coefficients s„ obtained from the fits are negligi- 
bly small and therefore only the coefficients Cn are shown 
in Fig. [TOl The major contributions to the AMR come 
from the terms with n = 2, 4, 6 and 8 indicating that the 
action of the Lorentz force on the mobile charges is not 
the source for the observed AMR in H-ZnO samples. 

As shown for the system ZnO-Cu with oxygen vacan- 
cies p^ . however, the influence of hydrogen cannot rule 
out the Zn-d contribution to the magnetic order as well as 
from the oxygen p-states. X-ray magnetic circular dichro- 
ism measurements are necessary to obtain the required 
information on the elements (and bands) contributions to 
the observed magnetic order. In particular the origin of 
the clear anomalous behavior at 50 if < T < 100 K with 
the unexpected change of sign of the AMR effect (factor 
C2, see Fig. 10) requires further studies that go beyond 
a magnetotransport characterization. We note, however, 
that a change of sign of the AMR (and the thcrmopower 
S") at T ~ 50 if has been observed for particular field di- 
rections in U3AS4 and U3P4 ferromagnetic single crystals 
[4^. The author interpreted the abrupt change and sign 
inversion of AMR (and 5") in the frame of spin-orbit cou- 
pHng (SOC) and a large sensitivity in energy of the spin 
density of states at the Fcrmi-level due to a spontaneous 
trigonal distortion in magnetically ordered state. 
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FIG. 9: Temperature dependence of the absolute value of the 
AMR defined in Eq. for the three H-ZnO samples. Note 
that actually the AMR changes sign at T ~ 50 K. 




50 100 150 200 

Temperature T (K) 

FIG. 10: Fitting coefficient Cn as a function of temperature 
obtained from the fits to Eq. ((9| to the experimental data 
shown in Fig. [S] The values of the coefficient C2 are divided 
by 4 to include them for clarity with the other coefficients. 
Note that it changes sign at ~ 50 K. 



F. Anomalous Hall Effect 

The anomalous Hall effect (AHE) was reported in 
magnetic- ion doped ZnO systems in the past j44l - |46l | but 
not yet in an (magnetic-ion) un-doped ZnO systems and 
at room temperature. The Hall resistance in ferromag- 
netic materials consists of two contributions which are 
the ordinary Hall resistance (due to Lorentz force) and 
the anomalous Hall resistance (due to an asymmetric 
scattering in the presence of magnetic order) and can 
be expressed by the following equation 



^Haii — Rh{H) + Rahe{M) , 



(10) 



where Rh and the Rahe are the ordinary and anoma- 
lous Hall resistances and M is the magnetization. The 
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FIG. 11: Anomalous Hall resistance of the H-2 sample as 
a function of a magnetic field at several temperatures. The 
linear background from the conventional Hall effect was sub- 
tracted from the measured curves. Note that AHE shows 
an anomalous change of field curvature above 50 K and it is 
nearly temperature independent above 225 K. 



dominant feature of the Hall data in our H-ZnO sam- 
ples is a linear dependence of i?Haii with magnetic field 
with a negative slope due to the ordinary contribution 
Rh- After subtracting Rh{H) from the measured data, 
an anomalous Hall effect contribution is obtained for all 
three H-ZnO samples, as shown in Fig. [TT] for sample 
H-2. Clear s— like loops with very weak hysteresis are 
observed in Rahe{H) up to 300 K. The shape of the 
loops at temperatures T < 50 K indicates that the hy- 
drogen related paramagnetic centers dominate, in agree- 
ment with SQUID results, see Fig.[lja). At intermediate 
temperatures 50 < T < 150 K the behavior of Rahe{H) 
is anomalous in the sense that it has a different field 
curvature without saturation at large fields. Note that 
in the same temperature range the carrier concentration, 
magnetoresistancc and anisotropic magnctorcsistance be- 
have also anomalously. At temperatures above 150 K the 
Rahe{H) curves follow the expected behavior for a fer- 
romagnet, see Fig. [TT] 

In order to further investigate whether the anoma- 
lous Hall effect is mainly affected by the magnetization 
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FIG. 12: A comparison of the SQUID and the anomalous Hall 
resistance data at 300 K. 



response of the samples we compared the SQUID and 
Rahe{H) results. Figure HU shows the M{H) and 
Rahe{H) loops at 300 K. Both curves are very similar, 
which indicate that the AHE originates from the intrinsic 
fcrromagnctism of the H-ZnO samples. 



crystals and in the available temperature and magnetic 
field range. The magnitude of the magnetoresistance in- 
creases with H-concentration. We observed the anoma- 
lous Hall effect and anisotropic magnetoresistance in the 
H-ZnO single crystals. The magnitude of anisotropic 
magnetoresistance was found to be 0.4 % at 250 K, a 
value comparable to polycrystalline cobalt. The anoma- 
lous Hall effect data showed a quantitative agreement 
to the SQUID data for all three H-ZnO single crystals, a 
fact that excludes impurities as the origin of the observed 
ferromagnetism. The observation of anisotropic magne- 
toresistance up to room temperature strongly suggests 
the presence of a spin-splitted band with a non-zero spin- 
orbit coupling in H-ZnO single crystals. At temperatures 
below 100 K, anomalous behaviors in the magnetoresis- 
tance, anisotropic magnetoresistance, and carrier density 
were observed that are apparently related to the change 
of the main contribution to the measured resistance, i.e 
from semiconducting to VRH-like transport. We believe 
that our findings would be useful for further understand- 
ing of defect induced magnetism in ZnO as well as other 
oxide systems and could be the starting point towards 
an efficient and reproducible way of inducing fcrromag- 
nctism in ZnO systems. 



IV. CONCLUSION 

We have studied the magnetic and magnetotransport 
properties of H- implanted ZnO single crystals with differ- 
ent hydrogen concentrations in the atomic percent range. 
Clear ferromagnetic-like loops were observed in all three 
H-ZnO samples at room temperature with a magnetiza- 
tion at saturation up to 4 emu/g. The Hall measure- 
ments confirmed the n-type transport mechanism in H- 
ZnO. We observed a negative magnetoresistance in all 
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